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Abstract. We review in details some recently proposed kinetic models of opinion dynamics. 
We discuss several variants including a generalised model. We provide mean field estimates 
for the critical points, which are numerically supported with reasonable accuracy. Using non- 
equilibrium relaxation techniques, we also investigate the nature of phase transitions observed 
in these models. We also study the nature of correlations as the critical points are approached. 



1. Introduction 

Application of statistical physics to understand social dynamics is an interesting and very active 
area of research [U [21 HI [5] . The key question one asks is how a set of interacting individuals 
choose between different options (vote, language, culture, opinions etc), leading to a state of 
'consensus' in one such option, or a state of coexistence of many of them. Opinion dynamics 
is one of the most important aspects of a society. It is a collective dynamical phenomena. 
Numerous models have been introduced so far to study the dynamics that leads to different 
opinion states and the processes that determine transitions between such states. Voter model 
which has a binary opinion variable with the opinion alignment proceeding by a random choice of 
neighbours [6] and the Sznajd discrete opinion formation model where more than just a pair of 
spins is associated with the decision making procedure [7j are two such examples. In some other 
models more than two opinions and also opinion as a continuous variable has been considered 
[HI El [TOl lllj . Models where the range of interactions are more than nearest neighbors but 
finite and even time dependent (see Ref |il2j and references therein) are also studied. 'Opinions' 
are subject to changes due to binary or group interactions, global feedback and even external 
factors. The usual interest in these studies lies in the distinct steady state properties: a phase 



characterized by individuals with widely different opinions and another phase with a measured 
fraction of individuals with similar opinions. 

We focus our attention to a specific class of models proposed recently |13tll4j. having apparent 
similarity with kinetic models of wealth exchange . The opinions of individuals are continuous 
variables in [—1, 1] which change due to binary interactions. The only parameter in these models 
is 'conviction', which is a measure of how much an individual sticks to his/her previous opinion 
while interacting with another. The system of such individuals, or the 'society', reaches a 
state of 'consensus' if this parameter stays above a threshold. One can also generalise this 
model [TB] by introducing another parameter modeling the 'influence' of the other individual. 
This study models the fact that the ability to influence need not be identical to one's conviction. 
When 'influence' and 'conviction' are identical, one gets the original model ^T3| I14j. The case 
of 'consensus' is an ordered phase (full or partial) while another 'disordered' phase exists, 
characterized by null (or very small) value of individual opinions. 

In this article, we will first review the original model and its generalisation. Then we will 
present some interesting features of the above model \13\ [H] , and propose some simplifications 
and variants. In particular, using Non-Equilibrium Relaxation (NER) [17|, we study the phase 
transition observed in these models. The associated critical exponents are obtained. Further, we 
propose a simpler variant of the original model to capture the minimum ingredients required to 
obtain the phase transition. Also the effect of global feedback is studied. The paper is organized 
as follows: In Sec. [2] and Sec. [3] we will discuss the original model and two of its variants. In 
Sec. m we discuss the general results, and the non-equilibrium relaxation in Sec. [5j Next, we 
discuss a generalized model in SecEl a map in Sec. [7] and also introduce a new model in Sec. [HI 
We conclude with a summary and discussions of our main results in Sec. [9l 

2. The LCCC model 

The basic idea of this model originated from a multi-agent statistical model of closed economy 
|15j where agents exchange a quantity w defined as wealth. Initially each agent begin with a 
certain amount of wealth Wi, i = 1,2, . . . , N such that Wi > and the total wealth W = Wi 
is conserved. The system evolves with a prescribed trading rule where agents interact with each 
other through a pairwise interaction characterised by a "saving" parameter A, with < A < 1. 
The dynamics (CC model) is as follows : 

w'i = Xwi + e(l - X){wi + Wj) 

w'j = Xwj + {l-e){l- X){wi + Wj) (1) 

where e is a stochastic variable that changes with time and as w is a conserved quantity, for each 
transaction w', + w', = Wi + Wj, where w', and w', are the agent wealth after the transaction. The 
functional form for steady state distribution fyw) is seen to be close to the T distribution |18] . 

Following this, Lallouache et. al. |131 [Ti] proposed a minimal multiagent model for the 
collective dynamics of opinion formation. Let Oi{t) G [— be the opinion of an individual i 
at time In a system of individuals, opinions change out of binary interactions: 

Oi{t + l) = \[oi{t) + eoj{t)] 

Oj{t + l) = \[oj{t) + e'oi{t)] (2) 

where e, e' are drawn randomly from uniform distributions in [0, 1] . Here, A is a parameter, which 
is interpreted as 'conviction'. The above model [131 E! (LCCC model hereafter) considers a 
society where everyone has the same value of conviction A. It is important to note that there 
are no conservation laws here. The opinions are bounded, i.e., —1 < Oj(t) < 1. The ordering 
in the system is measured by a quantity (order parameter) O = \ Oi\/N . Another important 



quantity is the so called 'condensation fraction' p, which is the fraction of the agent having 
Oi = ±1. Numerical simulations show that the multiagent system (dynamics given by Eqn. ([2])) 
goes into either of the two possible phases: for any A < Ac, Oj = Vi, while for A > Ac, O > and 
O ^ 1 as A ^ 1, with Ac ~ 2/3. Ac is the critical point of the phase transition. The relaxation 
time, defined as the time to reach a stationary value of O in time, diverges as r ~ |A — Ac|~^ 
when A ^ Ac on either side. A similar behavior is also observed for p. Although the values 
of the exponents differ [z ~ 1.0 it 0.1 and z ~ 0.7 it 0.1 with O and p respectively, reported 
in |14j). the critical points are same. The order parameter exponent /? is defined (for any order 
parameter X) as: 

X-(A-Ac)''. (3) 

For O its value is 0.10 it 0.01. Now if one also attempts to fit the growth of p in a similar form, 
one finds its value to be 0.95 it 0.02 . 

A mean field calculation can be proposed for the fixed point o*: 

o*[l-A(l + (6))]=0, (4) 

from which it is easy to show that the critical point is Ac = 1/(1 + (e)) ((. . .) refers to average). 
For uniform random distribution of e, (e) = 1/2 and hence, Ac = 2/3. It is important to note 
that this mean-field treatment does not incorporate the cut-offs at ±1 and yet gives the correct 
critical point. We also note that the underlying topology (Id, 2d or infinite range) has barely 
any effect on the critical point. 

It is appropriate to mention here that a map version of the model [ 131 E] has also been 
proposed: 

o(t + l) = A(l + e(t))o(t) (5) 

where e(t) is drawn randomly from a uniform distribution in [0, 1], and o{t) is bounded in [—1, +1] 
as before. The critical value Ac can be analytically shown to be exp[— (21n2 — 1)] ~ 0.6796 [14j. 

3. Two more models 

3.1. A simpler model 

A simpler model (model C hereafter) can be proposed: 

Oi{t + l) = \oi{t) + eoj{t) 

Oj{t + l) = Xoj{t) + eoi{t) (6) 

where e, e' are drawn randomly from an uniform distribution in [0, 1] as earlier. Here, A is a 
parameter as in LCCC model, which is interpreted as 'conviction'. An individual i upon meeting 
another individual j, retains his own opinion proportional to his conviction but also picks up a 
random infiuence of that of the other. Numerical simulation reveals that the system goes into 
either of two possible phases: for any A < Ac, Oj = Vi, a symmetric phase, while for A > Ac, 
O > and goes to 1 as A —t- 1, a symmetry broken phase, with Ac ~ 1/2. Ac is the critical point 
of the phase transition (Fig. [1]). 

A mean field calculation similar to the one discussed earlier, can be proposed for the fixed 
point 0*: 

o*(l-A-(6)) = 0, (7) 

from which it is easy to show that the critical point is Ac = 1 — (e). For uniform random 
distribution of e, (e) = 1/2 and hence, Ac = 1/2. 
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Figure 1. The phase diagrams for the 3 models. Left: Behavior of order parameter O. Right: 
Behavior of condensation fraction p. 

3.2. Global effect on LCCC model 

In the context of social opinion formation, global opinion often takes a vital role in influencing 
one's opinion. In that case, a person i, apart from being "influenced" stochastically by a person 
J, is also "influenced" stochastically by the average opinion of the entire society at that moment. 
Mathematically the dynamics can be represented by 

Oi{t + l) = \[o,{t) + eoj{t)] + e'0{t) 

Oj(t+l) = \[oj{t) + r]Oi{t)]+r]'0{t), (8) 

where e, e', r/ and r/' are random numbers, drawn from uniform distribution in [0, 1]. In this case 
(model G hereafter), the symmetry broken phase O 7^ appears for A > 1/3, and for A < 1/3 
the system is in a symmetric phase, with O, = Vi and all individual agents have the opinion 
(Fig.d]). 

This transition point can again be explained by a mean-field approach. At the steady state, 
i.e. when O reaches a steady value, Eqn. ^ can be written as, 

o* = A(l + (e))o* + (e')o* (9) 

from which it can be easily shown that Ac = 1/3. 

As in LCCC, the critical points of these two models have barely any effect due to change in 
the underlying geometry. In what follows we study the multi-agent dynamics of these models 
and their non-equilibrium relaxation. 

4. Results 

As mentioned above for the mean-field version of LCCC, the order parameter exponent is, 
0.10 lb 0.01. The order parameter exponent for the model C is found to be 0.17 it 0.01 for O and 
again if one fits the behavior of j* in a similar form it comes out to be 0.98 it 0.02, which is similar 
to the values found for LCCC. For the model G, we estimate the exponent /3 = 0.081 it 0.001 
and for p it is 0.85 ±0.01. 

We also calculate the relaxation behavior of the order parameters. An order parameter X 
relaxes to equilibrium in time t as 



X{t) oc exp(-t/r) 
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Figure 2. Behavior of relaxation time r with A — Ac from the order parameter O, for model 
LCCC (From Ref. [H). The system size was N = 500. 

for the symmetric phase (going to 0) and 



in the symmetry-broken phase, Xq being the equilibrium value of X. Away from the critical 
point, we plot this relaxation behavior for different values of the parameter A, to extract the 
value of r. We observe that the relaxation time diverges as r ~ |A — \c\~^ both below and 
above the critical point Ac, with roughly the same exponent z. For LCCC (mean field), our 
estimate for z is 0.97 =b 0.01 for O while it is 1.10 it 0.01 for p, which is consistent with another 
estimate (1.16 =b 0.03 in Ref. [16]) but is different from that reported in Ref. [H]. For model 
C, our estimates for z are 1.58 ± 0.01 for O and 1.34 it 0.01 for p. In case of model G, for O, 
z = l.2± 0.1 and for p, z = 1.75 ± 0.01. 

In the following subsection, we also report the study of a lattice version of the LCCC model. 
In this version, the agents are arranged on a Id lattice, and a randomly chosen nearest neighbor 
pair update their opinions according to Eqn. ([2]). The critical behavior of this model is studied 
in detail. 

5. Non-equilibrium relaxation 

NER is a well established simulation strategy to investigate the phase transitions of systems in 
the thermodynamic limit (for a review see Ref. \n\ IT9] ) . In this method the system is allowed to 
relax from an initial non-equilibrium state to the equilibrium state. It turns out that with much 
smaller amount of systematic errors, the critical exponents can be estimated from this method. 

All thermodynamic quantities show a relaxation in time, starting with initial values away from 
equilibrium. This temporal relaxation of the NER functions (above mentioned thermodynamic 
quantities) is used to determine the critical point accurately, as well as the critical exponents. 
The most important quantity in this study is of course the relaxation of the order parameter for 
the transition. We study the relaxation of both O and p at and near the critical point. 

Here we first study the lattice version of the LCCC model in details, later we get back to the 
infinite range model as well. The dynamics is started from an initial condition with full order 
(i.e., all agents have same extreme opinion, +1 or —1). Away from criticality, the system is 
expected to relax exponentially to its equilibrium order parameter value. At the critical point, 
however, the order parameter will relax asymptotically following a power law of the form 



X{t) oc Xo[l - exp(-t/r)] 



o{t) ~ t-^. 



(10) 
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Figure 3. The relaxations of the two quantities O (order parameter) and p (condensation 
fraction) are shown at the critical point (Ac = 0.66679 ± 0.00001) for the Id version of LCCC 
model. The power-law decays fit to exponent l.OOibO.05 for O and 1.15ib0.01 for p. The system 
size was N = 1200. 

Note that a similar power-law is also observed for the condensation fraction p. Fig. [3] shows 
the behavior of the two quantities at the critical point. From the slope of the log-log plot, the 
critical exponent 5 can be found. From this figure we observe that the critical point for both 
the quantities are same, and also the critical exponent 5 is very close for both these quantities 
{5 = 1.15 ± 0.01 for relaxation of p and 1.00 ± 0.05 for O). 

The critical point of the transition can be obtained with high accuracy by plotting the decay 
of the order parameter at and on either side of the critical point. Fig. S] shows the variation of 
the order parameter O for different A values near criticality. The accurately estimated critical 
point turns out to be Ac = 0.66679 ± 0.00001 for A'" = 1200, which is very close to the value 
(2/3) quoted in |141ll6j. A similar estimate was made with p (Fig. S]), the estimate of the critical 
point turns out to be exactly the same. 

From usual finite size scaling theory (1J7|, the order parameter is expected to follow a scaling 
relation of the form 

X{t)^t-^T{t^/''\\A), (11) 

where A = A — Ac and is a universal scaling function of a form such that for large argument, the 
time dependence drops out {T{x) ~ x^'^H). From the data collapse of the off-critical relaxation 
(Insets of Fig. 2]), one can obtain the correlation time exponent to be ~ 1.5 ± 0.1 for p and 
?a 1.2 ±0.1 for O. 

We carry out the above analysis for the fully connected version of the LCCC model. We only 
quote the main results of this version. The critical point is estimated to be Ac ~ 0.66659±0. 00002 
for = 1200, slightly different from the lattice version. The critical exponent 5 for p turns out 
to be 1.2 lb 0.1, again close to value with the lattice measurement. The off-critical scaling for p 
yields the correlation time exponent ~ 1.1 it 0.1, slightly different from the lattice counterpart. 

We repeat the analysis for the fully connected version of model C. For A'^ = 512, the critical 
point is found to be 0.500097 ± 0.000001, the same for the decay of O and p, close to the mean- 
field estimate (1/2). We also get estimates of decay exponent 6 as 0.500 ± 0.005 for O and 
0.521 ±0.005 for p. Using the above scaling relation for the data collapse (Eqn. fjlip . we estimate 
values of exponent as 0.10 ± 0.01 for O and 2.00 ± 0.02 for p (Fig. [5]). 

The same analysis is repeated for the fully connected version of model G. For = 512, the 
critical point is found to be 0.33338 ± 0.00001, the same for the decay of O and p, close to 
the mean-field estimate (1/3). The estimates of the decay exponent 6 are 0.585 ± 0.001 both 




Figure 4. Left: Relaxation of the order parameter O near Ac for the Id version of LCCC 
model. The parameter A differs by 0.00001 for the successive sets. The power law is obtained 
(central curve) for Ac = 0.66679. The system size for the simulation is = 1200. Inset: Plot of 
t{X — XcY^^ against 0{t)t^ using the previously obtained value of J. The data collapse is obtained 
for z^ll ~ 1.2ib0.1. Right: Same for the condensation fraction p. Inset: Plot of t(A — Ac)*^'! against 
p{t)t^ using the previously obtained value of 5. The data collapse is obtained for v\\ « 1.5 ± 0.1. 
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Figure 5. Relaxation of the condensation fraction p near Ac for model C. The system size 
for the simulation was = 512. Inset: Plot of t(A — Ac)*^" against p{t)t^ using the previously 
obtained value of 5. The data collapse is obtained for v\\ ~ 2.0 it 0.1. 



for O and p. From the above scaling relation for the data collapse (Eqn. (|ll|) . the estimates of 
exponent are 1.6 ± 0.1 for O and v\\ = 2.0 ±0.1 for p. It is to be noted that these exponents 
violate the scaling laws. 

In the quest of finding if there is a growing length scale in the system as the critical point 
is approached, we adopted the following strategy: we begin with all agents having opinion 0.0 
except a particular agent was fixed at the negative extremity (—1) and was kept fixed forever 
in all subsequent updates. Now as the system is allowed to relax, it is expected that the agent 
with the rigid opinion will have a neighborhood of influence. This is a measure of the correlation 
length = -^-/-^j where is the number of agents with opinion Oj < — e, where e is a very 
small number (in our simulation we have taken e = 0.005). This is expected to grow as the 
critical point is approached from the sub-critical regime, which is indeed seen from the Monte 



Table 1. Table comparing the different quantities for the 3 models 



Model 


Ac (Mean field) 


Measured quantity 
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0.97(1) 


1.00(5) 


1.2(1) 


P 


0.95(2) 


1.1(1) 


1.2(1) 


1.1(1) 


C 


1/2 


o 


0.17(1) 


1.58(1) 


0.500(5) 


0.10(1) 


p 


0.98(2) 


1.34(1) 


0.521(5) 


2.00(2) 


G 


1/3 


o 


0.081(1) 


1.2(1) 


0.585(1) 


1.6(1) 


p 


0.85(1) 


1.75(1) 


0.585(1) 


2.0(1) 



0.9 
0.8 
0.7 
0.6 

o 
jjlT 

0.5 
0.4 
0.3 
0.2 
0.1 

0.0001 0.0002 0.0003 0.0004 0.0005 0.0006 0.0007 0.0008 

Figure 6. The correlation as a function of (Ac — A) in lei LCCC model for N = 1000. 
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Figure 7. The schematic phase diagram of the generalised model. Distinct symmetric phase 
(SP) and symmetry broken phase (SBP) exist. The \ = ji line corresponds to the Model LCCC 
while ^ = 1 line corresponds to Model C. 

Carlo simulations and will remain one in the super-critical regime. In Fig. 6 we have shown the 
growth of length scale for the Id LCCC model. 

6. A generalised model 

A generalised version of the LCCC model was introduced by Sen in Ref. [16] where the 
influencing parameter of the iih. agent Aj is in general different from the conviction parameter 





Hi. Here the binary opinion exchange is mathematicahy represented by 



Oiit + 1) 
Oj{t + l) 



\iOi{t) + efijOj{t) 
XjOj{t) + e'niOi{t) 



(12) 



where the variables are as defined in previous sections. In the study A and h were assumed to 
be homogeneous. The special case of A = ^ is the LCCC model. The limiting case of // = 1 
corresponds to the model C (Fig. [7]). In the generalised case the transition from symmetric to 
symmetry broken phase is determined by both A and fi, the mean field phase boundary given 
by A = 1 - f,/2. 

6.1. Non-universal behavior 

Along the phase boundary mentioned above, the critical behavior is reported to be strongly 
non-universal. Here also the relaxation time diverges close to the transition points along the 
phase boundary as r ~ |A — \c\~^ for corresponding values of ^c- But in this case z varies with 
He rather systematically, indicating a non-universal behavior. For /ic = 0.4, z = 1.04 ± 0.01; 
He = 2/3, z = 1.10 ± 0.03; Hc = 0.9, z = 1.21 it 0.01. This non-universal behavior is also present 
in the order parameter exponent /3. For the order parameter O, /3 = 0.079 it 0.001 for Hc = 0.4 
and f3 = 0.155 ± 0.001 for Hc = 0.9. 

The condensation fraction p also shows similar behaviour as O. The relaxation time diverges 
close to the transition points along the phase boundary as Tp ~ |A — Xc\~^'' for corresponding 
values of Hc- The value of Zp now only weakly varies with fic but is very close to z, indicating 
the existence of only one time scale. Also p ~ (A — Xc)^'' with /3p « 0.91 for Hc = 0.4, ^ 0.95 for 
He = 2/3 and ~ 1.0 for Hc = 0.9. Here also, although the non-universality is present, it is very 
weak. 

7. A generalized map 

Recently, an interesting map has been proposed by Chakrabarti in Ref. [20], limiting cases of 
which correspond to known kinetic models. It has the following form: 



where qi, 02 and are linear functions of a single parameter a, (0 < CKj < 1 for i = 1,2,3), 
with —00 < n < 1 and 9 taking value either a positive, finite value (take 9 = 1, for convenience) 
or 00 (sufficiently large value) and et,£,t are uniform in [0, 1] and are independent. 

(a) When ai = a and 02 = = 1 — a, n = 1, 9 = 00, x(t + l) = (a + e(l — a))x(t) + ^(l — a) for 
which the steady state distributions of x are close to Gamma distributions for large alpha. 
For ^ = e, this reduces to x(t + 1) = ax(t) + e(l — a){x(t) + 1), closely representing the CC 
model. 

(b) When ai = a and 02 = 03 = 1 — a, —00 < 0, 6* = 00, then x{t + 1) = a + e(l — a)x{t) + 

n-2 

^(1 — a)". It is observed that x is distributed as a power law as P{x) = x "-1 . 

(c) When ai = 02 = a, = 0, n = 1 and 9 = 1, then Eqn. (fT3]l reduces to x{t + 1) = 
mina(l + e)x(t), 1, which is nothing but the map version of the LCCC model as given in 
Eqn. dS]). 

8. A model with preference 

Agents are endowed with opinion to begin with, with Oj distributed random uniformly in 
[— 1, +1]. The society is also parametrized by a conviction A (as in LCCC). 



x(t + l) 



min{(ai + eta2)x{t) + CtQ^, 9} 



(13) 
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Figure 8. Behavior or order parameter O for different system sizes. 

The interactions are binary. A pair of agents i and j meet and finds out who has the stronger 
opinion (compare values of Oj and oj). If Oj > Oj, then 

Oj = \{oj + eoi), (14) 

and if Oj < Oj 

Oi = \{oi + eoj), (15) 

e being the usual random number drawn from an uniform distribution in [0, 1] . Note that there 
is no dynamics if they have the same opinion (they agree). 

One computes the usual order parameter O = \ Oi\/N and also the fraction p of agents 
at ±1. O seems to undergo a phase transition at some Ac ^ 0.52 (Fig. [8]) while p is a step 
function at the same Ac- These signatures tempt us to conclude that the phase transition here 
is discontinuous in nature. 

9. Summary and Discussions 

We have reviewed here the studies on kinetic exchange model for opinion formation \13\ [H] , its 
two variants (proposed here) and a generalized model [TB]. In particular, we report the results 
from the study using non-equilibrium relaxation simulation technique to study the relaxation 
behavior of the order parameters for similar models. 

Here the lattice version of the LCCC model was proposed. The critical points and the critical 
exponents S and for Id version were found using two quantities, i.e., the average opinion O 
(order parameter) and the condensation fraction p. The critical points for both the quantities 
are exactly the same, however, the values of the exponents differ slightly. The mean field version 
gave same critical point and slightly different exponent values. Then a simpler version of the 
LCCC model was studied. And finally we also study the effect of the feedback of the global 
opinion formation in this model. A mean field calculation is made for both these variants, giving 
the critical points to be 1/2 and 1/3 respectively. These are again numerically verified within 
good accuracy. Analysing the NER functions of the order parameters, the dynamical critical 
exponent, the order parameter exponent were estimated. 

Next, we discuss the generalized model with 'influence' and 'conviction' |16j . and comment 
on the limiting cases which correspond to models LCCC and C. We also discuss the recently 
proposed map model, the special cases of which resemble some known kinetic exchange models 




in literature. The phase transitions observed in all the above mentioned cases has the unique 
feature that the disordered phase consists of individual opinions with values very near to 0. 
While in other phase transitions in the disordered phase the corresponding individual order 
parameters have finite positive or negative values but collectively resulting to a value of nearly 
0. We also propose a modified model in the same spirit as the models discussed earlier, the only 
difference being that the agents are only influential over the ones with weaker opinion. However, 
agents who agree, do not influence each other. The model shows clear characteristics of a first 
order transition. 
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